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SUMMARY 

Die three-dimensional nonlinear partial differential equation 
for the velocity potential in a supersonic stream is transformed hy 
the method of characteristics to obtain a system of three ordinary 
differential equations wherein all quantities are known or can "be 
calculated except the three velocity differentials. When the ordi- 
nary differential equations are converted to difference equations, 
the velocity components at an unknown point in a network of charac- 
teristics can "be calculated if the velocity components at three 
neighboring points of the network are known. Use of the difference 
equations for computing the supersonic potential flow past unsym- 
metrical boundaries is discussed. 

Application of the method of characteristics to the linearized 
three-dimensional equation results in a relatively simple system of 
difference equations that can be used to compute the supersonic flow 
past boundaries for which no other linearized solution is available. 


INTRODUCTION 

The velocity potential in a steady supersonic stream is described 
by a nonlinear partial differential equation of second order with 
three independent variables. If the velocity throughout the disturbed 
portion of the flow field is assumed to differ very little from the 
velocity in the undisturbed stream, the complete equation can be 
reduced to a linear equation whose form is identical to the wave 
equation of physics. Many solutions of the linearized equation have 
been obtained for particular types of boundary condition. The pres- 
sure distributions on thin wings, for example, can now be determined 1 
for many plan forms and profiles (references 1 to 5). Solutions are 
also available for determining the pressure distribution on axially 
symmetric bodies (references 6 and 7) or on cones that may be 
unsymmetri cal (reference 8). No linearized solutions are known, 
however, for general body forms or for general wing-body combinations. 
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When the Mich number and the velocity change appreciably as 
the flow passes a solid bounding surface, the complete differential 
equation must be solved if accurate results are required. For a 
few special cases, such as plane flow (reference 9) and axially 
symmetric conical flow (reference 10), the exact equation has been 
solved without use of the method of characteristics. These solutions 
were- obtainable because the problems to which they apply are com- 
pletely determined by a single independent variable. For plane flow, 
the local. Mach number is deter mine d by the angle of the streamline 
with respect to the base plane and, for axially symmetric conical 
flow, the local Ifech number depends only on the angle of the conical 
ray with respect to the axis of symmetry. When reduction of the 
problem to a single independent variable is impossible, solutions 
for the exact equation are available by application of the method 
of characteristics. 

In references 11 and 12, the method of characteristics is applied 
to the calculation of the supersonic flow in a field that is plane 
or has axial symmetry. For plane flow, the procedure is essentially 
a step-by-step development of the Prandtl -Meyer solution (reference 9), 
whereby the effect of interacting characteristics is determined. Each 
characteristic has a constant strength throughout the field and the 
inclination of the characteristic With respect to the local stream 
direction is the local Mach angle. For axially symmetric flows, the 
strength of a characteristic varies with the distance from the axis 
of symmetry; hence the problem cannot be reduced to a single vari- 
able. The computations are therefore somewhat more laborious than 
for plane flows, but the simplicity of geometry is maintained. The 
flow can still be completely represented in a single coordinate plane 
and the angle between the characteristic and the local stream direc- 
tion is still the local Mach angle. For plane and axially symmetric 
flows, the method of characteristics has been extended to include 
rotational flow (reference 12). 

When the flow variables vary simultaneously in all three 
coordinate directions, graphical representation is more difficult 
and numerical computations are more lengthy. The relation between 
the characteristic surfaces and the coordinate planes is no longer 
simple. A formulation of the procedure for unsymmetrical flows has 
been published by Ferrari (reference 13), who derived the ordinary 
differential equations that determine the variation of velocity 
components along intersections of characteristic surfaces with 
meridian planes and discussed the procedure for the case of an 
axially symmetric body at angle of attack. 
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A reformulation and generalization of the methods described in 
reference 13 was completed at the NACA Lewis laboratory in March 1948 
and is presented herein. Hie physical concepts used in the method 
of characteristics for general flows are described in some detail 
and the required difference equations are derived for the three com- 
mon orthogonal coordinate systems. With these equations and con- 
cepts, the potential flow past three-dimensional boundaries can be 
determined by point-to-point computation to any desired degree of 
accuracy. By application of this method to the linearized differen- 
tial equation for supersonic potential flow, solutions can be obtained 
for boundaries for which no other solutions are yet available or for 
which other solutions involve more computation than the procedure 
described. No attempt is made to determine the effect of shock waves 
on the flow field. All shock waves are assumed to be replaced by 
isentropic compression waves in order to retain the assumption of 


potential flow. To the extent that this assumption is valid, the 
resulting solution can satisfy the boundary conditions of the real 
flow field. 



SYMBOLS 

The following symbols are used in this report: 

A 


partial differential equation of characteristics 

Ai,A2,As 

Bi,B2,B 3 

C]l,C2,C3 

d 1> d 2> D 3 

Vs J 

> 

quantities depending on velocity components and 
coordinates at points of characteristic network 

A°,A* ,A". 

. .A m 

source of characteristic surface that extends 
downstream and away from solid boundary (type I) 

a i,k 


coefficients of D^^y in potential equation 

B°,B’ ,B" . 

. .BP 1 

source of characteristic surface that extends 
downstream and toward solid boundary (type II ) 

C°,C*,C". 

. .0“ 

source curve determined by intersection of two 
characteristic surfaces 

D i,k^ 

= 

^ 2 

(i,k = 1,2,3, table I) 

ox^ox^ 

D°,D» ,D". 

. .Dm 

source curve defined by intersection of charac- 


teristic of type II with = constant surface 
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E°,E' ,E" . 

■dm 

• ♦ 

source curve defined "by intersection of charac- 
teristic of type I with x^ = constant surface 


. .i® 

derived source curve on solid boundary and lying 
in X]_ = constant surface 

f(x x ,X 3 ) 


function defining integral surface 
S = x 2 - f (x^Xg) = 0 

<F,G' ,<?’ . 

. .cF 

intersection curve of characteristic of type H 
with solid boundary 



function in potential equation whose nature 
depends on coordinate system used 

k 


i 

ratio of critical speed to speed of sound 

M 


Mach number 

R 


ratio of coordinate differences required to solve 
simultaneous difference equations for velocity 
components 

r,e,q> 


spherical coordinates 

S 


integral surface of potential equation 

T 


dVjr dvj 

function of x ±) v ± , and 

1 3 

y 


ratio of local flow velocity to critical speed 

v i 


first partial derivatives of ¥ with respect 
to xj (i = 1,2,3) 

v x>V v z~l 

v x’ v r> v e 

v r' v e> v <pJ 

► 

ratio of velocity components to critical speed 
in Cartesian, cylindrical, and spherical 
coordinates, respectively 

x,r,cp 


cylindrical coordinates 

x,y,z 


Cartesian coordinates 



NACA TN No. 1849 


orthogonal coordinates used in derivation 
ratio of specific heats (7 = 1.40 for air) 


5 


Xi,X2>X3 

7 

6 

€ 


1 


X 

n 

¥ 

Subscripts : 


A°,A' ,A" . . .GF 


1,2,3 

Superscripts : 

m 

n 


angle between tangent to solid boundary and 
X 3 = constant surface 


angle between tangent to solid boundary and 
= constant surface 

factor required to eliminate coordinate factors 
from Vj[ to establish boundary condition 
(table I) 


dx? 

si 

c)Xo 

5^ 


along characteristic 


along characteristic 


ratio of velocity potential to critical speed 


values of quantity at point of characteristic 
network 

coordinate directions, x-^Xg^Xj 


last point of source curve 

general intermediate point of source curve 


THEORY 

Physical Concepts 

The physical basis of the method of characteristics is the 
nature of wave propagation in a supersonic stream. Consider a moving 
stream of compressible fluid that is uniform except for the effects 
of a single stationary point source of disturbance (fig. l). If the 
stream is moving subsonically, the wave front of a disturbance (shown 
at successive intervals of time in fig. l) has no stationary envelope j 
when steady-state conditions are obtained, the discontinuous wave 
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front has disappeared from the field of interest and a continuous 
disturbed flow field remains. If the stream is moving supersonically 
(fig. 1(h)), the wave front possesses a stationary conical envelope. 
The vertex of this envelope is the point source, its axis is the 
stream direction, and its conical half -angle is the Mach angle. 

When steady-state conditions are obtained, the initial spherical 
surface of discontinuity has disappeared, but the conical envelope 
of the disturbed portion of the field (Mach cone) remains. Die Ifach 
cone is a surface of infinitesimal discontinuity upon which the stream 
variables have two values, that of the free stream and that of the 
disturbed stream. Such a surface, upon which two solutions simulta- 
neously exist, is called a characteristic and its differential equa- 
tion can be directly obtained from the differential equation of the 
entire flow field. 

Consider now a space curve that acts as a source of disturbance 
in a supersonic stream (not necessarily uniform) (fig. 2). The shape 
of the space -curve source may be arbitrary except that each point 
of the curve is assumed to lie outside the Mach cones emanating from 
all other points of the curve. In figure 2 the curve is assumed, 
for convenience, to lie in a plane perpendicular to a coordinate axis 
and the flow is in the general direction of the x^-axis. The 
envelope of all disturbances from this source forms two distinct sur- 
faces. These surfaces are characteristics and the stream variables 
may have either of two solutions at any point on the surfaces. If 
the line source is regarded as a series of point disturbances, the 
shapes of the characteristics are clearly such that they form an 
envelope to the Mach cones emanating from all point sources on the 
curve. If only very small regions of the flow field are considered, 
the characteristics are closely approximated by their tangent planes. 
These tangent planes and their intersections with coordinate planes 
and with each other are used in the method of characteristics to con- 
struct continuously varying flow fields (reference 13). (For spher- 
ical coordinates, the surfaces X]_ = constant become spheres rather 
/ than planes, but because most problems are best adapted to the use 
of Cartesian or cylindrical coordinates, the term "coordinate plane" 
is sometimes used to describe a surface for which either x^ or 
X3 Is constant.) 

The initial sources of the characteristics are actual or 
assumed discontinuities at the boundaries of the flow. These bound- 
aries consist of the body past which, or through which, the flow is 
to be computed and the surface that separates the disturbed and 
undisturbed regions of the flow field. Smoothly varying boundaries 
are replaced by a succession of tangent planes whose spatial dimen- 
sions are small so that the original boundary is closely approximated 
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■by the resulting polyhedral surface. For boundary sources, it is 
convenient in most cases to choose curves that lie in surfaces for 
which one of the coordinates is constant. If Cartesian or cylindri- 
cal coordinates are used, for example, all initial sources can he 
so chosen that they lie in Xj = constant planes. If the flow vari- 
ables are known on an initial source lying on the solid boundary, a 
characteristic surface of type I (fig. 2) is determined. Intersection 
of this surface with a characteristic of type II, whose source is a 
«n».n distance from the initial source, determines another curve that 
nay be regarded as a secondary source of two more characteristic sur- 
faces (one of each type). Hie secondary source so determined does 
not, in general, lie in an = constant plane, but a procedure is 
subsequently described whereby all secondary sources can be transferred 
to planes parallel to the initial source. 

When each characteristic is replaced by its tangent planes, the 
flow field is represented by a network of points, each of which is 
a Junction of two or three lines of intersection of planes tangent 
to the characteristic with coordinate planes. Hie mathematical theory 
of characteristics provides equations that determine the velocity 
components at these Junction points when conditions at neighboring 
Junction points are known. 


Mathematical Development 

Hie theory of characteristics is presented in many texts dealing 
with differential equations (for example, references 14 to 16). Hie 
following development considers the theory for nonlinear partial dif- 
ferential equations of the second order in three independent variables. 
Hiis type represents the equation for compressible potential flow in 
a three-dimensional orthogonal coordinate system. Although the fol- 
lowing development differs from the derivation given for cylindrical 
coordinates in reference 13, the final total -differential equations 
obtained are identical. 

Hie differential equation of the velocity potential for com- 
pressible, nonvlscous, irrotational flow can be written as 

a li D n V + a 12 I) 12' if + a 13 D 13 W + a 22 D 22 T + a 23 D 23 V + a 33 D 33 ¥ + H = 0 

( 1 ) 

where a i ^j c and H are functions of velocity components and coor- 
dinates and the symbols Di # fc¥ represent the second partial deriva- 
tives of ¥ with respect to the coordinates. Hie nature of the 
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quantities In equation (l) for Cartesian, cylindrical, and spherical 
coordinates is indicated in table I. The irrotational -flow condi- 
tions required to derive equation (l) are 


■>12’*' " 

B 13f ' »31* ' (2) 

“as' 1 ' ■ V 

Hie feasibility of integrating equation (l) depends on the 
existence of surfaces in the flow field along which the derivatives 
of W are continuous. If discontinuities are possible in the field, 
these surfaces (often called integral surfaces) must be so defined 
that the flow variables may be discontinuous across the surfaces but 
not along them. The differential equation that defines such a family 
of surfaces is derived as follows: Assume that the surface 
S(x^,X 2 ,x 3 ) = 0 is an integral surface of equation (1). If S is 
solved for one of its independent variables, such as 13 , then . 


S = ig - f(x 1 ,x 3 ) s 0 (3) 

and the relation between the Increments dx^,dx 2 > and d x 3 on 
this surface is 


where 


dxg — XdX]_ + ndx 3 



(4) 


dx? 

--SR 


The differential of v 1 may be written 


4T i ' D 1,1 “*1 + + B i,3 fdI 3 


(i = 1,2,3) (5) 
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If one of the coordinate Increments , for example dx 2 , is eliminated 
from equation (5) "by means of equation (4), then equation (5) becomes 


dv l s < D i,l T+ XD i,2 y ) ta l + 


( D i,3^ + •*®l,2^ dx 3 


( 6 ) 


Because x^ and are Independent variables, the increments 

dxj and dxj may be alternately set equal to zero to obtain the 
relations 


*1,* - fe) - M, i,2 ? (7) 

' b constant 




constant 


- >*1,** 


(7a) 


where 
x 3 = 


dv4 dvj 

t— - and are taken 

^3 

constant planes and x^ = 


along the Intersections of S with 
constant planes, respectively. 


With equations (7) and (7a), all second derivatives of W 
except one can be eliminated from equation (l) j thus, if D 22 F is 
the exception, equation (l) becomes 


AD 22 ¥ + T = 0 


( 8 ) 


where T is a function of Vj, x^, 
coefficient of D 22 W is 




and the 


_ 2 2 
A ■ a^X - + a^fy + a 22 - a 23 n + aj^i (9) 

Equation (9) is called the characteristic form, of equation (1) . If 
A ^ 0, then D 22 ^ and consequently all other partial derivatives 
of V are uniquely determined by equation (8) and are consequently 
single -valued. If A = 0, however, the partial derivatives are not 
uniquely determined by equation (8) and may therefore be multivalued. 
Hence the possibility that discontinuities may exist in the flow 
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field is established by the condition A = 0. Die integral surfaces 
defined by the relation A = 0 are called the characteristics of 
equation (l). If the characteristics exist and are real, they con- 
stitute a family of surfaces along which equation (l) is necessarily 
integrable. Other surfaces may cross the characteristics and hence 
contain discontinuities in the derivatives of ¥ . In general, if the 
flow is supersonic, real characteristic surfaces exist. An interest- 
ing exception is the equation for conical flow, which is discussed 
in appendix A. 

If the coefficient of (equation (9)) is set equal to 

zero and the resulting equation is solved for A , the differential 
equation of the characteristics of equation (l) becomes 

- ( a 12 ” 1^13) ± V (^ a 13 “ a 12^ 2 " ( a 22 “ MS -23 + ^ 2& 33 ) 

x ^ 

(10) 

where the + and - are used for characteristics of types I and 
H, respectively (fig. 2). Die quantities A and n may be called 
intersection parameters of the characteristic surfaces because they 
determine the rate of change of X 2 with respect to xq and X3, 
respectively, along the characteristics. (For plane and axially 
symmetric flower the velocity vector can always be represented in a 
single plane. Die parameter u is then zero and equation (10 ) can 
be reduced to 

A = tan ( 0 ± p) 

where 0 is the angle between the local velocity vector and the 
x-axis and p is the local Mach angle.) 

For three-dimensional flow the significance of A and n is 
shown in figure 3, where x^, xg, and X 3 are represented as 
Cartesian coordinates. If in figure 3(a) the intersection of the 
characteristic with an x^ = constant surface is given (line A°A' ), 
then two values of A at point A° can be calculated from equa- 
tion (10), in which the approximate value of p at A is given 
by the difference ratio ^ 

= ( ~ X 2,A' ~*2,A° 

c0n8timt ‘ ^.-3, AO 
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These two values of X determine the slope of the intersection at 
A 0 of a characteristic of each type with the x 3 = constant plane 
containing A°. If the intersectio n of a characteristic with an 
13 » constant plane is given (line AOC°, fig. 3(b)), two values of 
U can be calculated from the equation obtained by solving the charac- 
teristic equation (10) for |i . Because the x^-axis is assumed to 
be in the general direction of the flow, the parameter u will 
usually be given and the two values of X will be computed from 
equation (10). 


The expressions for the changes in the velocity components in 
passing from one point to a neighboring point along a characteristic 
intersection can be determined by setting the function T, obtained 
in equation (8), equal to zero. Because the nature of this function 
has not yet been specified, the derivation will be repeated in mor e 
detail. With the aid of the irrotational flow conditions (equa- 
tions (2)), the following relations are obtained from equations (7) 
and (7a): 


D n ¥ 

B i# 

B 2# 

B 3# 


dvi dv? „ 

Gf - * asf ♦ ’'’V' 


dvc 


JET - V 


dv-a 

Si 


- * D 23 ¥ 


dv, 

_ ^22^ 


dv 


dv 9 


dx* “ M dxT + 


22 




> 


and equation (8) becomes 


( 11 ) 


T 


" a ll dv l + ( a 12 - Xa ll) dv 2 + 

[ dv 2 dv, 

(®23 " * a 13 " ^33) dx“ + a 33 + H 

3 3 


dXjL = 0 


( 12 ) 
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where the coefficient of Dj > 2 V has been set equal to zero in accord- 
ance with the characteristic condition. Equation (12) provides a 
relation between the velocity increments in moving a small distance 
dx^ along intersections of characteristics with x 3 = constant 
surfaces. If, for example, the velocity components and coordinates 
are known at A° and A' in figure 3(a), all quantities in equa- 
tion (12) are known except dvn, dv 2 , and dv 3 . With the aid of 
equations (7) end (7a), the relation 

dv 2 dv 2 dv 3 dv x 
A dx 3 * * dx 2 + ax 2 " dx 3 

is obtained and equation ( 12 ) can be converted to 


a 13 dv l + (*23 ” |ia 33 )dv 2 + a 33 dv 3 


av 2 av^ 

+ (a 12 - Xa!! - ua 13 ) 5 ^ + a n ^ + H dx 3 = 0 


which relates the velocity components in moving a distance dx 3 along 
Intersections of characteristics with x^ = constant planes. For 
cylindrical coordinates, equations (12) and (13) are identical to 
those derived by Ferrari (reference 13). When equations (12) and 
(13) are converted to difference equations, the supersonic flow past 
arbitrary boundaries can be constructed by a point-to-point process 
of numerical Integration. 

Formulation of difference equations along characteristic inter- 
sections. - Bie ^problem of constructing the supersonic flow past 
arbitrary boundary surfaces may be divided into two parts: (l) Given 

two neighboring source curves, both of which lie in surfaces for 
which Xi = constant, determine a third source curve lying in another 
x^ = constant surface slightly downstream; ( 2 ) given one source 
curve on the solid boundary and another curve just off -the solid 
boundary, find another source curve (sli^itly downstream) on the 
solid boundary. (All source curves are again to lie in surfaces for 
which xq = constant.) These problems and the geometry of the 
solution are illustrated in figure 4, where x-j_, x 2 , and x 3 are 

represented as Cartesian coordinates. 

Suppose that the points A°, A', A", F°, F’, F", G°, G', and G" 
are located on a body that represents a boundary of the flow, ihe 
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rectangles AOA'F'F° and A' A'^'F* are assumed to be sufficiently 
awftii to closely approximate the actual contour of the body in the 
region A^F'F 0 . Die rectangles A^'E'E 0 and A'A'^'E* are 
tangent planes that closely approximate a characteristic surface of 
type I whose source is the curve A°A' A" on the body; the rectangles 
BOB'D'DO and B’B"D"D’ are tangent planes that closely approximate 
a characteristic of type II whose source is B°B' B" . It is desired 
first to determine the coordinates and velocity components at points 
E°, E' , and E" (or D°, D' , and D" ) when the coordinates and velocity 
components are known at points A 0 , A’ , A" , B°, B’ , and B" . When 
•the new field source E°E'E" has been determined, the second problem 
requires the determination of the velocity components at the new 
surface source F°F' F " . 


Bie fact that lines such as A'E 1 and B'D' do n ot ge nerall y 
intersect at the same value of Xj_ as the lines AOE° and B°D° 
mast be noted in solving the first problem. If the body is very 
unsymmetrl cal , the xi value of the intersection point changes 
appreciably as the successive points E' ,E" . . .E® are computed. 

Such lines as A’E' and E°E' , however, which are intersections of 
the same characteristic tangent plane with x 3 = constant surfaces 
and xi = constant surfaces, respectively, naturally intersect at 
the chosen values of x^ and x 3 . 

The velocity components along the derived source curve E°E'E" 
can be determined with equations (12) and (13), which may be written 
as follows: Along an intersection of a characteristic with an 
X 3 = constant plane, 

a ll dv l + ( a l2 ” ^an)dV 2 + a 13 dv 3 + P^dx^ = 0 (14) 


where 


dv 2 dv 3 

P 1 = ( a 23 * Xa 13 ' ^33^ dx^ + a 33 dx^ + H 


Along an intersection of a characteristic with an x^ = constant 
plane, 

/ 

a 13 iv l + ( a 23 " ^33^ dv 2 + a 33 av 3 + ^3^3 = 0 


(15) 
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where 

dv 2 dv-^ 

p 3 = ( a 12 " (^13 " Xa ll) + a ll + H 

If the velocity gradients and the coordinates change only 
slightly in moving from one point to a neighboring point on the 
network of characteristic intersections, the differentials in equa- 
tions (14) and (15) may be replaced by differences. The coordinate 
differences dx^ and dx 3 are easily obtained. For any point such 
as C° (or C’ or C" ) , which is a junction of two intersections 
of characteristics with an x 3 = constant plane, the two relations 


* X 2,C° ^,B° 

A B° = 5T 2 — _ x ■ ■ 
1,C° 1,B° 


and 

v 

A Xi, CO - X lfAO 

can be used to determine x^ qO and x 2 go (x 3 qO is given for 
such points). Explicitly the required coordinates are 


X 1,C° 


X 2,C° 


X 2,A° - X 2,B° + X B° x l,B° ~ X A° x l,A° 
X B° “ X A° 

\° X 2,A° ~ ^A^.B 0 +X b° X A° ( X 1,B° 
X B° X A° 



(16) 


For such points as E' , D’ , E" , and D" , which are junctions of 
intersections of the same characteristic tangent plane with an 
x-[_ = constant plane and with an x 3 = constant plane, the coor- 
dinates Xjl and x 3 are known from preceding points and x 2 can 
be determined from relations such as, for point E' , 


- x 2,E« " *2,^ 
A ' *1,E« " X 1,A’ 


( 17 ) 
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In order to determine the velocity components at points such as 
E' f three equations are required. From, equations (14) and (15) the 
following difference equations may he obtained: 

Along A’E' 


( a ll) A , < V * 1,E' “ V 1,A' ) + ( a 12 " Xa ll^ A , ( V 2,E» - V 2,A' ) 


+ ^ a 13^ A , ^ V 3,E* ' V 3,A' * + P 1,A« ^E' " X 1,A' ^ = 0 

(18) 

Along B’D* 

^ a ll^B» ^ V 1,C’ “ V 1,B') + ( a 12 - ^ a ll^B* ( V 2,C' " V 2,B*) 

+ ^ a l3^ B , ( V 3,C’ " t 3,B‘) + P 1,B' ( X 1,C' " X 1,B» ) = 0 

Along E°E» 

( a 13^E° ( V 1,E' “ V 1,E°) + ( a 23 " ^33^10 ( V 2,E* “ V 2,E°) 

+ ^ a 33^jo ^ V 3,E' “ V 3,E°) + P 3,E0 ( X 3,E* " X 3,E°) “ 0 

( 20 ) 

Because only two characteristic intersections Join at points 
such as E' , three independent relations between the three velocity 
components at E' cannot be obtained unless the assumption that the 
quantities v A vary linearly for small distances along characteristic 
intersections is used" to eliminate v ± c . from equation (19). This 

assumption is expressed by the relation 

V 1,C’ ~ V 1,A' = *1,0' ~ *1,A’ = 

Y i,E’ " V i,A' X 1,E’ “ ^A’ 


E 


( 21 ) 
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With equation (21), equation (19) "becomes 

( a ll^ B . fr^E 1 " V 1,B* + (l-B) v l,A»] 

♦ (a 12 - l>2,E- - T 2,B' * (1 - K)t 2,A'J 

+ ^ a 13^-g, C Rv 3,E* " V 3,B' + ( 1 " E ^ v 3,A t ] 

+ *1,B' ( x 1,C* “ X 1,B’ ) ° 0 ( 22 ) 

Bie solutions of equations (18), (22), and (20) are 



(23) 


where 



A = Ai (B 2C3 

- b 3 c 2 ) - 

a 2 ( B 1 C 3 " B 3 C 1^ + A 3 

, ( b iC 2 - 

B 2 Ci) 

*1 

= (a U>*. 

B 1 

= (a 12 - ^ a n ) A , 

C i 

H 

II 

h. 

- ( a ll) B ' 

®2 

= (a 12 - Aa n ) B , 

C 2 

= ^ a 13^ B » 

A 3 

% 

0 

j? 

n 

B 3 

= (a 23 - ua^)^ 

C 3 

= ^33)^ 


D 1 = ” A l v l, 

A* ~ B 1 v 2, 

,A' ” C 1 V 3,A' + P 1,A' 

< X 1,E' - 

X 1,A' ) 
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®2 - TT (*2 T 1 

,A« H 

h B 2 v 2,A* + 

1 

C 2 V 3,A* 

J 

) 


- 5 ( A 2 T 1 

,B' H 

h B 2 V 2,B' + 

\ 

C 2 V 3,B' 

) + R (x l,C' 

** 

d 3 = " a 3 v 1,E° " 

B 3 V 2,E° " C 3 V 3,E° 

+ P3,e°( x 3,E^ " 

Z 3 

[ 

- 

dv* 


dVo 1 


P 1,A' = 1 

- a 33 

d^ + (a 23 

- Xa 13 ■ 

-^33) d^ + =J A , 


__ 

dv. 


dv 2 


P 1,B' = 

_ a 33 

5^ + ( a 23 

- Xa 13 

" ^ 33 ) d^ + H - 

B' 


r 

dvi 


> dv 2 

1 

P 3,E° = ' 

L a ii 

dx” + (a 12 

- na 13 

- a ia>sr + H . 

E° 


The difference forms of the derivatives in the expressions for 
P-, »i, Pt , and Pt_ go are given hy equations (B16) and (B17) 
in , appendix B. ’ 


By means of equation (23), the velocity components at any point 
I ? 1 of a derived source curve can he determined provided that the 
velocity components at points A n , B * 1 and E&-1 are known. Inas- 
much as the source curves A°A' < . .A m and B°B' . . .B® are given, 
the velocity components at all A and B points can he considered 
known. Ihe determination of the E° point, however, is also required 
before the computation can proceed. Because only two equations are 
available to determine the velocity components at E°, (one along 
APE° and one along B°E°) , one of the velocity components at E° 
must he known or assumed. If the problem has a plane of symmetry, 
that plane can he used as a reference plane on which the cross com- 
ponent of the velocity V 3 is zero. For most practical problems, 
such a reference plane can usually he found. If none is available, 
however, the computations must proceed by trial and error, that is, 
an initial value must he assumed for one of the velocity components 

at E° and the new source E°E' . . .IP 1 must then he computed entirely 
around the body. If the final computed value of the velocity com- 
ponent is considerably different from the assumed value, this value 
must he correspondingly modified and the computations repeated. Bie 
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resulting process nay be much too lengthy in practice. Die value 
of V3 in the reference plane Xj = 0 is therefore assumed to he 
known and equal to zero. ETam equations (14) and (15) , the following 
difference equations for the velocity components at E° are then 
obtained : 

( a ll) A o < V 1,E° - V 1,A°> + ( a 12 * Xa ll> A o < V 2,E° ’ V 2,A°) 

+ P;l,ao ( X 1,E° " x l,Ao) 55 0 ( 24 ) 

( a n) B ° (vi,E° " V 1,B°) + ( a 12 ' Xa llV ^ t 2,EO - V 2,A°) 

+ Px,B° (x 1#E o - Xi,a°) " 0 ( 25 ) 

whose solutions are 

A 1 D 1 
d 2 

a 2 b 2 

where 

Ai = (a U ) A o ' B i - (<^ 2 - SlV 

h. 88 ( a ll)so B 2 = ( a 12 - ^llh 0 

D 1 = “ A 1 v 1,A° “ B 1 v 2,A° + P 1,A° ( X 1,E° ■ X 1,A°J 
D 2 = - A 2 v-^ b o - BgVg^o + Px^ B o (xx^jo - Xx^ A °) 


®1 D 1 

B 2 D 2 
Vb! 

a 2 b 2 


(26) 
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f dv 2 dv 3 

1,A<? “ 

\ "^33 dx 3 + a 33 dx 3 


( dv 2 dv 3 

1,B° = 

\ ^33 dx 3 + a 33 dx 3 


B c 


This solution can also be obtained from equation (23) by setting 

A,, = 3J= D » 0 and C = R = 1. 

3 3 3. 3 


With equations (26). and (23), the new field source 

E°E’ . . .E® can be computed. Biis new source curve determines 
two new characteristics, one of which extends toward the body 
whereas the other extends away from the body. Hie characteristic 
that extends away from the body may Intersect another character- 
istic of type II whose source is a curve above B°B’ . . .TP- 
in figure 4. Bale intersection can be calculated by the procedure 
used to calculate E°E' . . .EP; thus the entire characteristic 
whose initial source is A°A' . . .A m can be confuted if con- 
ditions on a surface a small distance upstream of it are known. 

Ihe determination of a new Initial source on the body, such as 
F°F* . . .1®, must now be considered. 


The X]_ -coordinate of the surface in which the new initial 
source F°F' . . .I® 1 is to be located may be chosen as the 
Junction of the type-H-characteristic intersection from E° 
with the body. Ihls value of Xj F o is determined from the 
relations ’ 


tan 5po 


1 

^2,A° ( X 2,FQ “ x 2,Ao) 
Z 1,F° " X 1,A° 


X - ,E° 

E° = x X ,pO - x ljE o 


where tan 8po is the slope of 
equation (10) with 


A°F° and 


X 2,E’ " ^.E 0 

Ueo = =-* — 

X 3,E° “ X 3,E° 


Xpo is calculated from 
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Die factor l/t) 2 Is given in table I for each of the coordinate 
systems. 


The coordinates of points that are junctions of the character- 
istic intersection from E' with the body can be determined from the 
relations 


tan 8 Qt 





_ X 2,A' ) 
Z 1,A» 




X 2,G» ~ X 2,E« 
X 1,G» ” X 1,E' 


For determining the velocity compone nts at F* , two char acter - 
istic equations are available, one along E*G' and one along F°F* . 
Die third independent relation is obtained from the condition that 
the velocity vector at any point on the body must be in a plane tan- 
gent to the body at that point. Diis condition is evidently given by 
the relation (fig. 5) 


T l2 v 2 = r l3' v ‘3 co ^ € + "Vi tan 5 (27) 

where e is the angle between the x 2 -direction and the trace of 
the body tangent plane in an xi = constant surface and 5 is the 
angle between the xi -direction and the trace of the tangent plane 
in an X3 = constant surface . Die quantities rjg and t|j are 
required to eliminate the coordinate factor from v* when angular 
coord inates are used (table I). Along the characteristic source 
F°F’ , the following relation is obtained from equation (15) : 


(ai3) p0 (▼ljF* ~ V 1,F°) + ( a 23 “ ^33)50 ( V 2,F’ " V 2,F°) 


+ (’ a 33^pO ( V 3,F* " V 3,F°) + P 3,F° ( X 3,F’ 


X 3,F°) s 0 

(28) 


Equations (27) and (28) are two of the relations required to determine 
the velocity components at F* . In order to obtain the third relation, 
an assumption s imilar to that made- in determining E 1 is required 
to e limina te vi^g* from the characteristic relation along E’ G* . 

Die required assumption is 
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▼1, G' “ V 1,A' _ X 1,G' ~ X 1,A' = R 
v i,F' " v i,A’ X 1,F’ " X 1,A’ 

When equation (29) is used to replace in the relation 

( a ll^E« ( V 1,G’ " V 1,E’^ + ^ a 12 ~ Xa ll^E' ^ T 2,G' ‘ V 2,E') 

+ ^ a 13^E' ^ V 3,G’ ‘ V 3,E' ) + P 1,E' ( X 1,G' " X 1,E* ) = 


the third required equation becomes 


( a ll) E , [^ljF* _ V 1,E* + (1 -H)vi,a*] 

+ (a 12 - ^ a 11 ) RI [e^F* ‘ y 2,E' + v 2 

+ ^ a 13^E* [^3^* “ V 3,E' + ( 1 "E)v'3^a»] 



(32) 

where 

A = A 1 ( b 2 C 3 “ b 3 C 2^ " a 2 (®1®3 " B 3 c l) + A 3 (®l c 2 ” B^l) 

A 1 = ' fcan 5 F' B 1 = - T l2,F* ^1 = (^3 co '*' c )pi 

a 2 » ( a 13) F o B 2 = ^ a 23 ~ ^ a 33^F° ° 2 = ^ a 33^pO 

a 3 = tell^E' B 3 = ( a 12 " Xa ll^E» C 3 " (*13^* 
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D 2 = “ A 2 v 1,F° " B 2 v 2,F° " C 2 V 3,F° + P 3,FO ( X 3,F‘ “ X 3,F°) 
D 3 = IP ( A 3 v 1,A‘ + B 3 V 2,A* + C 3 V 3,A') 

" £ ( A 3 V 1,E» + B 3 V 2,E' + C 3 V 3,E') + ^ ( X 1,G' “ X 1,E' ) 

[ dv 2 cLv-l I 

f a !2 ' Ma 13 - Xa ll) d^ + a ll d5^ + ^ 

P 1,E» = [(*23 * Xa 13 _ ^ 33 ) d4 + a 33 + ^ 

For point F°, which lies in the reference plane, v 3 is 
assumed to "be zero and equation (27) "becomes ° 


(^2 V 2^F° = V 1,F° ^ 5 F° ( 33 ) 

Die second relation required to determine v^po and v^po is 
obtained from equation (14), which, for V3 = 0, "becomes 


where 


^ a ll^E° ^ V 1#F° “ v l,Eo) + ( a 12 " Xa ll^E° ( y 2,F° " V 2,E°) 

+ P 1,E° ( X 1,F° " X 1,E°) = 0 (3*) 


? 1,E° = [ a 33 (as| “ * fcf) + H ]eO 


Die solution of equations (33) and (34) for v-^po is 

tell^go V 1,E° + ( a 12 " Xa ll^jo V 2,E° ’ P 1,E° ( X 1,F° “ X 1,E°) 
V 1,F° “ 


( a ll >]50 + (a 12 “ Xa ll^EO f “- 2 - 


po 


( 35 ) 
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With equations (32), (33), and (35), new initial sources on the 
surface can he computed when conditions on a field source and on 
another initial source (both slightly upstream) are known. 

Equations (23), (26), (32), (33), and (35) constitute the 
solution to the two basic problems mentioned at the beginning of 
the section. With these equations, the entire disturbed flow field 
about a body placed in a supersonic stream can be calculated pro- 
vided that (l) the flow remains everywhere supersonic, ( 2 ) the 
rotation and the viscosity of the flow are zero or negligible, 

(3) conditions are known along the surface that separates the uni- 
form flow field from the disturbed flow field, and (4) a reference 
surface exists for which V 3 = 0 . Condition (4) is a practical 
limitation that will eliminate trial -and -error procedure. It may 
be replaced by some other known condition satisfied by one of the 
velocity components at points of types E° and F° if the equations 
for these points are appropriately altered. 

The initial surface required by condition (3) may be, for example, 
the shock surface attached to the nose of a pointed body or to the 
lip of an open-nosed body. For unsymmetrical bodies, the determina- 
tion of the form and intensity of the initial shock surface is diffi- 
cult and the replacement of the shock by a characteristic surface 
may be essential for solution of some problems. This procedure 
neglects not only the rotationality of the real flow field, but also 
the abrupt compression at the foremost boundary of the disturbed 
portion of the field. The resulting solution therefore fails to 
satisfy exactly the real boundary conditions. ' For such problems, 
the linearized characteristic equations (to be derived subsequently) 
may yield a solution as accurate as that obtainable with the non- 
linearized equations. 

If the flow over a pointed body is to be computed, some assump- 
tion must initially be made concerning the velocity components near 
the point of the body because the vertex itself is singular. Perhaps 
the simplest assumption is that the flow is conical for a smai 1 dis- 
tance beyond the vertex. Any reasonable assumption is valid, how- 
ever, because the effect of that initial assumption can be made neg- 
ligible by scaling up the body so that the portion near the vertex 
is a very small part of the entire body. The initiation of the char- 
acteristic solution for pointed bodies is discussed in more detail 
in appendix A. For open-nosed bodies the starting process is some- 
what more straightforward because the velocity components a mrwll 
distance from the lip are accurately obtainable from two-dimensional 
flow theory. 
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The equations required to carry out computations are summarized 
in appendix B, whereas in appendix C the first steps involved in 
the computation of the flow past a cone that has an elliptic cross 
section are described. 

Solution of the linearized equation by means of characteristic 
surfaces. - When the changes in the velocity components throughout 
the disturbed flow field are small relative to the critical velocity, 
equation (l) can be reduced to a linear equation, which, for 
Cartesian and cylindrical coordinates,' has the form 

aiiDll^ + a22D22^ + a33D33V + H = 0 (36) 

where 

a n = (1-#) 

and where the quantities k^v^ 2 are neglected in the expressions 
for & 22 > a 33> and H (table I). The value of M is now assumed 
to remain constant throughout the flow field. 

The linearized theory for supersonic flow consists in finding 
solutions to equation (36) for various boundary conditions. For 
some types of solid boundary, such as thin wings and conical bodies, 
a variety of solutions are available. These solutions permit the 
determination of flow parameters at the surface of the boundary with- 
out determining the flow in the entire disturbed field. Solutions 
that have been obtained for bodies whose contour varies in the 
stream direction or for body-wing combinations postulate some 
restriction as to the shape of the boundary and in themselves involve 
considerable computation. For bodies that are not axially sym- 
metric, solutions have as yet been obtained only for the condition 
that the flow is conical (reference 8). The following discussion 
develops a method of characteristics for solution of the linearized 
differential equation for general boundaries . The resulting equa- 
tions are quite simple relative to those obtained for nonlinear flow, 
although the method involves, as with nonlinear flow, confutation 
of the entire flow field. The method is consequently not intended 
to be used when boundary-surface solutions are available unless these 
solutions involve more confutation than the one presented herein. 

The differential equation for the characteristic surfaces of 
equation (36) is 
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which, for Cartesian coordinates, 156000163 


(38) 


In equation (37), X and p have the same significance as in equa- 
tion (4); hence equations (7) and (7a) may he used to eliminate the 
partial derivatives from equation (36) in the manner used for non- 
linear flow. Die resulting ordinary differential equation that relates 
the velocity components in moving from one point to a neighboring 
point along an intersection of a characteristic with an x 3 = constant 
surface is 


= ± J 1 + 

yM 2 - 1 




dx^ + Hdx^ 


0 


(39) 


Die corresponding relation along intersections of a characteristic 
with an xi = constant surface becomes 


/dvi 


a 33 ( dv 3 - mAt 2 ) + a ll (dxj 


_ X 



+ Hdx 3 = 0 


(40) 


Equations (39) and (40) can also he obtained by setting a 12 , a 13 , 

and a 23 equal to zero in equations (14) and (15). Die difference 
equations for each type of source point may consequently be derived 
simply by setting these coefficients equal to zero in the nonlin- ' 
earl zed equations. Die resulting equations are given in appendix B 
for a general field-source point E& (equation (B 6 )) and for a 
general solid -boundary source point I ® 1 (equation (B7)). If v 3 
is again assumed to be zero in the X 3 = 0 surface, the velocity 
components for points E° can be obtained from equation (B 6 ) by 
setting D 3 = ^ E n-1 = 0 and B = 1. Die velocity components for 
point F° are given by equations (B 8 ) and (B9). Die expressions 
for the coordinates of the various types of source point are -the 
same as for nonlinear theory and are given by equations (B12) to (B15). 


With the expressions given in appendix B, "the linearized 
solution for the supersonic flow past general boundaries can be 
obtained provided that, to the conditions mentioned for the nonlinear 
solutions, the additional condition is added that the velocity is 
nowhere greatly different from the velocity of the undisturbed 
stream. Die initial surface that is required to start the computations 
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la, for linearized theory, the foremost characteristic surface of 
•Uie disturbed field, upon which it may be assumed that the velocity 
components have their free -stream values. An initial source curve 
an the boundary surface must be assumed, but the effect of this 
initial assumption can again be made negligible by scaling up the 
initial section of the body. 


SUMMARY OF ANALYSIS 

Die method of characteristics has been applied to the non- 
linear and linearized partial differential equations for the velocity 
potential in a supersonic stream." By use of the resulting difference 
equations, the velocity components can be computed throughout an 
irrotational supersonic flow field for arbitrary boundary conditions. 
Die solution for the linearized equation involves considerably less 
computation than the solution for the nonlinear! zed equation, although 
both require the determination of the entire flow field that influ- 
ences the flow variables at the boundaries. Die use of the linearized 
solution is suggested in those problems for which the disturbances 
resulting from the presence of boundary surfaces are everywhere small 
and for which no simpler linearized solution is available. The non- 
linearized solution is required, for accurate results, when the 
presence of solid boundaries results in large changes in the velocity 
components . No procedure has been given for treating the effects 
of an initial shock surface; however, the replacement of the initial 
shock surface with a characteristic surface should not lead to serious 
error unless the shock is intense. 


Lewis Flight Propulsion laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 29, 1948. 
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APPENDIX A 

TREATMENT OF UNSYMMETRICAL CONICAL FLOW FIELDS 

If the supersonic flow over pointed bodies is to he determined, 
the initial portion of the body near the point can usually be 
regarded as a small cone. If the nose is axially symmetric, the 
velocity components at the surface of the initial cone and through- 
out the conical portion of the field can be calculated by the method 
of Taylor and Maccoll (reference 10) . If the initial cone is not 
axially symmetric, however (if it has, for example, an elliptic 
cross section), no exact solution other than the method of character 
istics is known for determining the flow throughout the conical 
field. If a linearized solution is desired, the procedure presented 
in reference 8 may be used for determining the velocity components 
at the surface of an unsymmetrical cone. This method may be used 
to obtain the initial source curve required to start the solution 
by means of characteristics if the conical nose is 'only a small part 
of the entire body. If the conical nose is a considerable part of 
the body, however, flow conditions must be accurately known over an 
entire surface ahead of the characteristic whose source is the first 
variation of the body from the conical shape of the nose. A method 
for determining the velocity components throughout an unsymmetrical 
conical flow field must therefore be developed. 

The method of characteristics cannot be simplified to solve 
this problem because the assumption that the flow is conical leads 
to an equation that has no real characteristic solutions. If the 
conical -flow conditions are imposed on equation (l), then for spher- 
ical coordinates, equation (l) becomes 


a 22 D 22^ + a 23 D 23^ + a 33 I) 33^ + - 0 (A 1 ) 


where 


H‘ = (2 - k^g 2 - k^ 2 ) + cot 6 (1 + k 2 ^ 2 ) 

The characteristic form of equation (Al) is 


A = a 2 2 


H a 23 


^ 2& 33 


(A2) 
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If equation (A2) is set equal to zero, the characteristics can "be 
shown to he imaginary unless 


Because Vg and v^ are generally quite small near the cone surface 
the inequality is invalid throughout a conical field. 

The following procedure may he used, however, to determine the 
flow past unsymmetrical cones hy means of characteristics. Hie con- 
ditions at an initial source curve near the point can he fixed at 
some reasonable values, possibly based on a linear variation of 
velocity components with the conical half -angle 0. Hie normal pro- 
cedure described in the text can then be used to develop the flow 
field until there is no further variation of the velocity components 
along a given radius vector from the vertex of the cone. Hie number 
of computations required to establish the conical solution will 
depend on the accuracy of the initial assumption. Hie use of this 
procedure to compute the flow over an elliptic cone is described in 
appendix C. Other procedures, such as the direct numerical integra- 
tion of equation (Al) in holograph planes or along prescribed conical 
surfaces, appears to be feasible and their development should con- 
stitute a profitable field for future research. 



NACA TN No. 1Q49 


29 


APPENDIX B 


SUMMARY OF EQUATIONS 
Nonlinear Theory 

The following equations are used for field-source points E 11 : 


r i,nP ” ~ A 


T z,& 


T Z,'E a 


(Bl) 


(B2) 


A l( B 2^" B 5 B 2)" A 2( B 1^5~ B 5 B l)' tA 5( B l B 2" B 2 B l) 


-A 


f A 1 V 1 > E n+B l v 2 ,E n + °1 


(B3) 


where 


A - A^ (B 2 C3-B3C 2 )-A2(B 1 C3-B3C 1 )+A3(B 1 C 2 -B 2 C 1 ) 

A 1 = ^*11^11 ^ " ^ a 12" Xa U^ A n 

a 2 “ ( a ll)tf* *2 = ^ a 12" A ^1^ 

A 3 = ^ a 13^E n "^ % = ^ a 23”^ 


C 1 “ ^ a 13^A n 

C 2 = ( a 13) B n 
C 3 = (‘^e 11 " 1 


_A l V l,A n -B l V 2,A n ' C l V 3,A n + ^'l,A n ^ X l,E n _X l,A n ^ 
/ \ 

1-R 


D 2 = ~ VWlfA 10 - + B 2 V 2,A n + C 2 v 3,A n 


l.B 11 


’ 1 (Vl + B 2 V 2 + C 2 V 3 ) ^n + B (x l,C n “ x l,B a ) 
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= - ( A 3 V 1 ■* 

■ % v 2 + ( 

' 

- 

dv_ 

o 

P l,A n = 

a 33 


'■ 

- 

dV 3 


a 33 

d^3 


r 

dv. 

? 3,E n-1 = 

a n 


dv. 


a 33 dXj + ( a 23" A a 13 _|ia 33^ dx 3 + 


H 


»n 


dv, 


B 


n 


i ) ^ + = 


E 31 " 1 


x l t C n ~ x l,A n 
' X l,E n ' X l,A n 

For point E° (v,. = 0) : A^ = Bj = D 3 = 0, C 3 = E = 1. 


Equations for the velocity components v^ j,n, Vg F n, and 
Vj j, n at the solid-boundary source point F n have the same forms 

as equations (Bl), (B2), and (B3), respectively, where the coef- 
ficients are given by the following expressions: 


A i = tan V 


*2 = (* 13 ^ 11-1 
A 3 = ^ a ll^E n 


D 1 = 0 


®1 = - T l2,F n 

= ( a 23“^ a 33^n-l 
B 3 = ( a 12" X a ll^ 


C x = (tj 3 cot e) Fn 
C 2 = ( a 33^ F n-l 
C 3 = ^ a 13^E n 


D 2 = - (Vl + B 2 v 2 + C 2 v 3y j>n-l + \ F n-l < x 3 ,F n “ 

^5 * TT (*SV + ^3 V 2,A n + C 3%Aij 

- f (*3 T 1 + % t 2 + C 3 v 3) E n + (x l,G n ' X l,E n) 
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3,pn-l " 

dv 2 dv^ 

(ai 2 -U a l 3 -X a ll) ^ + an + H 

F n-1 

*1,*“ = 

dvg dv 3 

dxj + a 33 + H 

E n 


p X l»C? l ” X ltA 11 
X l^' X l,A n 


i F° (v_=0), the following equations are used: 

O 

( a ll) E o v ljE o+(ai2" Xa ll ) E o V 2 1 E 0 ' P l,E o(X l,F 0 ' X l 


1,F° 


( a ll) E o + ( a 12^ a ll) E o ( tQ 5J 2 & ) FO 


(B4) 


1,E° 


/ d ^3 

k 3 


dv<- 


a 33ydx-7 dx. 


+ H 


E^ 


^2*2^0 = ( V 1 tan 6 ) f o 


(B5) 


Linearized Theory 

The folloving equations are used for field- source points eF 1 : 

? 

^pip2~ ^gn^l 

v i,E n = ^(y^) 


V D i 


? 2,E n “ A^X^TX^y 


> (B6) 


^n-l^" 1 ^ 

= ‘ ^JVV" = ^ + V- 1 *2,1“ 
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where 
A x = M 2 -! 

®l = • A l' V l- XT 2 ) A n + P l,A' 1< \,E I ‘' X l,A n) 

1 P An pD 

Do = ^ An(v n -A _.v ) - r=- (Vn-XVp) + — ^ — (x -x ) 

2 E 1 l,A n B 11 2,A n E 1 2 B n E l,C n l,^ 1 


D 3 - (v 3- |1 ' , 2 ) E n-l * P 3 (E n-l (X 3,E“' Z 3,E n - l) 


l,A n 


**33 


/ dv 3 

V*3 ' 11 


dv 2 V 


+ H 


JA 


n 




a 33 


3 dv 2^ 

v 3 ^ " ^ 


H 


B* 


sX" 1 


X ^\ _H_ 

a 33 W ^7 a 33 


’V- 1 


The following equations are used for boundary- surface source 
points Jp: 


l,F n 


2,^ 


*1^(13 cot e \£% [\, f d~V- 1 (,| S cot C) pn 


A 3 D 2 (ti 3 cot e)^ n “I >3 tan 6 ^ 


>(B7) 


Al^C^g F n'^n tan 5 F n) _ ^n-l % tan ^n 

V 3,F» * 

(^^n-Ypn taD V 
(n 3 cot e) Fn 


J 
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where 


4 ‘ A 1 <T| 3 COt £) P» 'V ^ V] 

A 1 = M 2 -! 


D 2 " (vs-txVg)^! + 


A 1 / dv l ^ dv 2\ H 
*33 l toj_ to i/ a 33 


jOl-l <X 3,F n X 3,f J1 - 1> 


“3 " ¥ A 1 ( 7 l,An-V T 2 ,AD) - r < v l- Xv 2> E n + ^ 



dv- dv c 


P 1 " a 33 lds 3 " * dsjy 


+ H 


Je 11 


For point F° (v 3 =0), the following equations apply: 


(v 1 -Xv 2 ) 


^ (x 


E° A 1 1 


,F°’ X 1,E°^ 


1,F° 


1 -X 


( tan s\ 

E°\ ^2 A 


(B8) 


1,E° 


/ dv 3 dv 2 \ 


(VsU - (T i tan S) po 


(B9) 
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Expressions Involving Coordinates and Derivations of 
Velocity Components (Nonlinear or Linear Theory) 
The difference equations for p and X are: 

x ,-x ''N 

= 2 ,A n+1 2,A n 

A * X 3,A n+1 “ X 3,A n 
-x 


= 


2.B n + 1 2.B 11 


^ X 3,B“ +1 ^jB 11 


X - 

= 2.E 11 2,E n ~ J ~ 

^ 1 X 3,EP’ X 3,E n " 1 


(BIO) 






Nonlinear, 


/ 

-(ai 2 -P a i 3 )± y(pa 13 -a 12 ) 2 -4aii( a 2 2 -pa 23 -Hi 2a 3 3) 


A = 


-2a' 


11 


Linear, 


XBll) 


X = ± 


' a 22 + ^ a 33 


- a U 


where 


v 2 , V z 2 , V r 2 , Vg 2 and v^ 2 are neglected in a i#k 


J 


For computing E 11 , + is used for X n and X^n-i and - for 




For computing F 11 , + is used for X^.q and - for X^ n . 

The expressions for the coordinates of the points in the net- 
work are as follows: 
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Coordinates of C 11 are 


_ _ X 2.B n ~^.A nt V Z l.A n ' X B Ii:I l.B 11 

1 , 0 ” VV 


A 


X = x _+X (x -x ) 
2,C n 2,A n A n l,C n l,A n 


> (B12) 


X „ = X _ = X „ 

5,d°- 3,A n 3,B°- 


J 


Coordinates of E 11 are 


*l,E n = X 1,E° X 1,C° 




X 2,E n Z 2,A n+X A n ^ X l,E 11 " X l,A n ) 


> (B13) 


X = X = X 

3^ 3,A n 3,^ 


J 


Coordinates of G n are 

x r>-x „+x _ tan 6 -Ax 

x _ 2,E^ 2,A l,A n G 11 E n l t E n 

\ )G n = tan S^-X^ 


A 


x = x +A (x -x ) 
2,(? 1 2,E n E n 1^ l,E n 


} (B14) 


x = x 
3,G n 3,E n 


J 
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Coordinates of F n are 


x = x = x 
1,1® 1,F° 1,G° 




= x +(tu) n tan 6 (x -x ) ) (BIS) 

>,F n 2,A n 2 A n F 31 1,1® l,A n ( 


l 3,F n = 


J 


The velocity derivatives in difference form are as follovs: 
For P, , 


/ dV 5\ _ V 3,A n+1 ~ V 3 t A 11 

’ V +1 ‘V 


'dw 


2,A n+1 2 ,A n 






> (B16) 


J 


Similar expressions apply for points B 11 and E 11 . 


For P^, 


dv- 

L dx« 


. V l,E n ~ 1 ~ V l,A n ~ 1 

T / X -X 

1 V - 1 1, I ® -1 - * n " 1 


'dv. 


[\ Jjd 


1>A 


n-l“ v n-1 

2, A 


1 * n-l”* n-1 

x E n_1 1>E • 1,A 

Similar expressions apply for points F 11-1 




> (B17) 


J 
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APPENDIX C 


COMPUTATION OF FLOW OVER CONE 
WITH ELLIPTIC CROSS SECTION 

As an example of the use of the nonlinearized method of char- 
acteristics for unsymmetrical bodies, the first steps in deter- 
mining the flow over an elliptic cone are described in detail. This 
example illustrates the starting procedure for pointed bodies in 
general, because the shape of the body can be modified progressively 
downstream once the initial conical field at the point has been 
determined. 

The elliptic cone used is similar to that described in refer- 
ence 8 and is shown in figure 6. The ratio of the major to the 
minor axis of the cross section is 3 and the slope of the cone in 
the plane cp = 90° is 0.315; thus the cone is represented by the 
equation 


| = tan 6 = — (Cl) 

f\fs cos*tp+sin 2 cp 

and the expression* for cot e is 

cot e = — = 8 sin cp cos cp 

r< ^ 9 cos 2 cp + sin 2 cp 

The free-stream Mach number was assumed to be 1.90, 
ence 8. 

A preliminary computation plan is illustrated in figure 7, 
where each point in the plane cp = 0 (fig. 7(a)) represents a 
source curve composed of 17 points circumferentially distributed 
at values of cp of 0, 10, 20, 30, 40, 50, 60, 65, 70, 73, 76, 79, 
82, 85, 87, 89, and 90° (fig. 7(b)). As yet no basis other than 
experience is known for determining the number and distribution 
of points required for a prescribed degree of accuracy. It is 
evident, however, that a denser distribution of points is required 
where the shape of the surface varies most rapidly. 

In order to start the computations, values of v x , v r , and 

vcp were assumed at the points of source curves 1 to 6 (figs. 6 
and 7). At group 1, on the cone surface, the velocity components 
were assigned according to the equations 


(C2) 

as in refer- 
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y = v tan 8 + v cot e_ 

r,n x,n n cp,n n 


V^=v ^+v ^+v ^ 

v n v x,n + v r,n + v cp,n 


(C3) 


wherein V and v^ were assumed at each of the 17 points and v r 
and v x were then calculated. For n=0 (cp=0°), V was assumed 

to have the value resulting from two-dimensional compression of 
6 q = tan"^ 0.105; for n = 16 (cp = 90°), V was assumed to have 
the value obtained on a symmetrical cone of half-angle 8 lg = 

tan"^ 0.315. For n = 0 and n = 16, v,p must be zero because 
the flow is symmetrical with respect to these planes. For other 
values of n, the values of V and Vq> were selected by drawing 
reasonable-looking curves between their values at qp = 0° and 
cp = 90°. 

At the Mach cone (group 6), v cp = v r = 0 and v x = 1.585 (for 

the assumed free- stream Mach number 1.90). For groups 2 to 5 the 
velocity components end the values of r /x were assumed to vary 
linearly (for each value of n) between their values at the cone 
surface and at the Mach cone. 

With the velocity components and coordinates thus assigned 
at groups 1 to 6, the velocity components and coordinates at 
groups 7, 9, 13, 19, and 27 can be simultaneously computed by 
equations (Bl) to (B3) and (B10) to (B13) in appendix B. In order 
to compute group 7, for example, the points in group 1 are used as 

the A n points and the points in group 2 are the B n points . When 

group 7 (E°, E^, . . . E* 6 ) has been computed, these points 
together with group 1 can be used to compute the new solid-boundary 
source (group 8). 

An alternate procedure to the one outlined would consist in 
assigning values only at groups 1 and 6 and computing groups 31 
and 36. This procedure has the advantage that values need be 
assumed only at one surface source curve and at the Mach cone, 
conditions at the Mach cone being known if no shock is assumed to 
occur. The increments of r between the groups, however, are 
very large relative to the value of r at the solid boundary; 
thus the assumption that the coordinate increments are small is 
violated. It is therefore doubtful whether the computed values of 
the velocity components would converge as rapidly to their conical- 
field values as with the procedure previously described. The 
alternate procedure has the additional disadvantage that only one 
source curve can be computed at a time. 
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Because v^ is zero "both for cp = 0 and for cp = 90°, either 

of these planes can be used as the reference plane (fig. 7(b)). 

When two or more reference planes are available, any cumulative 
error involved in the computation of a source curve can be reduced 
by using each reference plane as a locus of E° and F° points 

and computing values of E 11 and F 11 from two directions. This 
procedure was used in the present preliminary computations. The 
difference between computed values of the velocity components 
resulting from the use of this procedure rather than computation 
in one direction is shown for group 9 in figure 8. Part of the 
difference between the two results is due to the difference in the 
value of x for the two groups of points, inasmuch as the n = 0 
points determine the x- coordinate of the entire group. Most of 
the difference, however, probably results from the fact that the 
values of some of the total derivatives, when approximated by their 
difference form, depend somewhat on the direction of computation 
(equations (BIO) and (B16)). 

In order to determine whether the number of points used in a 
group greatly affects the computation results, group 7 was com- 
puted twice. For the first computation 17 points were used and 
for the second 31 points. In addition, the assumed values of 
in groups 1 and 2 were slightly altered to determine the 

sensitivity of the computed results to changes in initial values. 
The results of these computations are shown in figure 9. The 
assumed values for groups 1 and 2 are indicated with solid lines. 
The alteration in the assumed values of v C p for the 17-point 
computation is indicated for group 1 by a broken line. The 
oscillation of the computed value of v x about the assumed 

initial values occurred for both sets of computations, and no 
difference in results is apparent that could be attributed to 
the use of a different number of points per group. The abrupt 
changes in v x at 0° and 20° for the 31-point group were attrib- 
utable to computation errors. These errors were retained to 
determine their effect on succeeding points in the group. The 
results indicate that the influence of such errors is limited 
to the immediate vicinity of'' the point for which the error is 
made. 


The difference in the initial values of v for the two 

cp 

groups resulted in relatively large changes in the computed values 
of v^ and smaller changes in the computed values of v x and 

v y . In the regions for which the assumed values of v^ were the 
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same (the leftward- running computations), the computed results were 
almost the same except near the peaks of the oscillations in v x 
and Vcp . ' 

The cause of the oscillations in the computed values of v x - 
and Vjp has not been definitely determined. The fact that an 

increase in the number of points did not reduce their amplitude 
indicates that the magnitude of the coordinate differences was not 
responsible. This conclusion is further substantiated by the fact 
that the oscillations occurred regardless of the direction of com- 
putation, although the peaks were somewhat displaced. The hypoth- 
esis is therefore advanced that the initial values chosen for the 
velocity components were unsatisfactory and tended to initiate 
strong compression and expansion waves in the flow field. Insuf- 
ficient time was available to check this hypothesis by altering 
the initial values. If the hypothesis is correct, a more accurate 
knowledge of initial values for each of the velocity components 
would appear to be very desirable for the computation of unsym- 
metrical flow fields by the method presented. Convergence of the 
computations may be very slow if the assumed values can be cor- 
rected only by means of strong compression or expansion of the 
flow. 


With ordinary calculating machines, a group of 17 points com- 
prising a source curve could be computed and checked in about 
80 man-hours. The computation setup consisted of 186 col umns for 

the E 11 groups and 130 columns for the F 11 groups. For a 
linearized solution these figures may be reduced by about one 
half. The use of electronic calculators will, of course, reduce 
both computation time and the possibility of error. 
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(a) Subsonic stream. (b) Supersonic stream 

Figure 1. «■- Illustration of propagation of initial disturbance in 
stream moving at subsonic and at supersonic velocities. 



Figure 2. - Illustration of two possible characteristic surfaces 
originating at source curve. 
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tan ^ X. 


tan X 


AM M 


V 


ban ^A(+) 

-tan-1 ^ A( _) 


(a) Velocity components 

given at points A 0 and A*. 


(b) Velocity components 

given at points A° and C°. 


Figure 3. - Illustration of method of determining characteristic inter- 
sections when velocity components at two points are known. 


In 

ll 

"pi 


C°, D°, E° 












Figure 4. - Network of characteristic intersections . near a solid boundary. 
All surfaces replaced by successive tangent planes. 
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ti 2 v 2 = ti^Vj cot € + v 1 tan 6 


Figure 5. - Illustration of boundary condition on velocity components 
Plane tangent to solid boundary at point P. 


-Initial assumed 
source curves 


wmmm 


Figure 6. - Sketch of elliptic cone used in example (appendix C), 
Free-stream Mach number, 1.90; a/b = 3.0; a/x = 0.315. 
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(a) Plane <p = 0, showing source curves assumed and to be 

computed. 

n = 10 9 8 6 420 

1 2 3 4 5 6 7 8 10 12 14 16 

1 1 1 1 1 1 — I — I. 1 1 ll 1 I 6 (Mach cone) 




1 (Solid cone) 





20 40 60 80 100 

<p t deg 

(b) Plane x = 1.0, showing distribution of points in 
assumed source curves. 


Figure 7. - Illustration of preliminary computation plan. 
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<P# de 8 

Figure 8* - Comparison of two methods of computing group 9 from 
assumed values in groups 2 and 3* 
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Figure 9. - Computed velocity components for group 7 from assumed 

values in groups 1 and 2. 



